Superconducting Proximity Effect through a Magnetic Domain Wall 
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We study the superconducting proximity effect in a superconductor-ferromagnet-superconductor 
(SFS) iieterostructure, containing a domain wall in the ferromagnetic region. For the ferromagnet 
we assume an alloy with an exchange splitting of the conduction bands comparable to the super- 
conducting gaps. We calculate the modification of the density of states in the center of the domain 
wall as a result of the proximity effect. We show that the density of states is sensitive to do- 
main wall parameters due to triplet-pairing correlations created in vicinity of the domain wall. We 
present a theoretical tool which in a very effective way enables retaining the full spatially dependent 
spin-space structure of the problem. 



PACS numbers: 74.45.-|-c,74.50.-|-r 

Introduction: Most promising candidates for meso- 
scopic devices with novel functionality are hybrid struc- 
tures containing superconducting elements. The key phe- 
nomenon that controls the behavior of such systems is 
the proximity effect. When a superconducting material is 
placed in contact with a normal metal (N), the supercon- 
ducting pair correlations leak over to the normal-metal 
side, changing its conduction properties in the vicinity 
of the separating interface. Quite similarly, the prop- 
erties on the superconducting side are also changed (the 
energy gap Aq is suppressed) due to the contact to a nor- 
mal metal. An alternative but equivalent way of thinking 
about the proximity effect is through Andrcev-reflectioni 
processes: an incoming electron from the normal side 
is transmitted together with another one as a Cooper 
pair into the superconducting side. This phase-coherent 
electron-hole conversion results in a nonzero pair ampli- 
tude in the normal metal. 

In the diffusive limit, the correlations relating to an 
incident electron with an energy E (the range of energies 
being set by the temperature T) above the chemical po- 
tential extend a characteristic distance of £,n = \/D/E 
into the normal metal)^ here D is the diffusion constant 
in N. If the extent of the N region is finite, another en- 
ergy scale, Et ^ D/L'^, enters the problem; L denotes 
the width of N. This so-called Thouless energy has as- 
sociated with it one of the generic features of diffusive 
superconductor-normal metal heterostructures, the mini- 
gap: the density of states in the normal metal develops 
a gap around the chemical potential in a manner similar 
to a superconductor (S) but with a smaller magnitude. 

If the normal conductor is replaced by a fer- 
romagnet (F), a multitude of new effects arise 
due to the emergence of yet another energy scale, 
that of the exchange splitting J of the two spin 
bands. Both on the theoretical2ilS*SiL2i24i2iii and on 
the experimental;^^i^^i^"^i^^i^^i^'''i^^i^^i^° side interest has 
grown recently in the rich physics of such systems. One 
source for new behavior is that, in the case with a singlet 
superconductor, the induced pair amplitude in the ferro- 
magnet is oscillatory.^^ However, the exchange splitting 



also gives rise to dephasing which, in turn, results in the 
decay of induced correlations over a characteristic dis- 
tance ^,7 = y^D/{E + J)w^ Unfortunately, since J is of 
the order of the Fermi energy Ep 'm typical ferromagnetic 
metals, this distance is very short. Still, experimental in- 
dication of the oscillatory behavior has been obtained 
in thin ferromagnetic layers and, relevant to the present 
paper, in weakly ferromagnetic alloys with J -C Ep^^^ 
Another question of current interest in SF proximity 
systems is the role of equal-spin triplet correlations 
If created, e.g. near magnetic inhomogeneities, such cor- 
relations would not be affected by the exchange splitting 
but could penetrate considerably longer distances into 
F(S2i Finally, the importance of domain walls has also 
been stressed for the Andreev conductancei2SiS&. 

In this paper we study an SFS structure, shown 
schematically in Fig. ^ in equilibrium. The ferromag- 
netic region consists of two domains with magnetizations 
oriented in opposite directions. The domains are sepa- 
rated by a domain wall, where the magnetization rotates 
continuously between the asymptotic values. While vary- 
ing in direction, the magnitude J is assumed constant 
throughout the F region. We show, that the local den- 
sity of states (LDOS) in the F region is strongly modified 
by the presence of the domain wall. In particular, we 
show that it can be very sensitive to the thickness of the 
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FIG. 1: SFS structure with two magnetic domains oriented 
along the z axis and separated by a domain wall of widt h dw ; 
dp denotes the length of the F region and ^ = ^ D/Aq is the 
superconducting coherence length. 
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domain wall in a certain parameter region. 

Basic equations: Proximity effect is a spatially inho- 
mogeneous phenomenon. An appropriate theoretical tool 
to treat such a problem is the quasiclassical theory of 
superconductivity,^^'^* which in its diffusive version has 
been formulated by Usadel^In equilibrium, the physical 
information is contained in the retarded Green functions 
G{z,E). Here, we assume spatial dependence in the co- 
ordinate z only, and E denotes the energy as measured 
from the chemical potential. The 4x4 matrix structure, 
arising from particle-hole and spin degrees of freedom, is 
denoted by the hat {") accent. 
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The off-diagonal elements determine the superconduct- 
ing pair amplitude. Quantities denoted with the "tilde" 
are related to those without one through A{z, E) = 
A{z,-E*)*. All the elements in Eq. (QJ are 2x2 spin 
matrices: e.g. Q = Qap with a,(3 = The Green 

fimctions satisfy the Usadel equation, 
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where the symbol (» denotes matrix multiplication, and 
-B] = A® B - B ® A- In writing Eq. |5J we have fol- 
lowed the standard way to describe ferromagnetic mate- 
rials through a spin-dependent energy shift which has 
the form Et^ — > Et^ — J ■ a. Here, T3 denotes the third 
Pauli-matrix in Nambu space, the vector J denotes the 
effective exchange field of the ferromagnet, and A is the 
superconducting order parameter (appropriate for weak- 
coupling spin-singlet pairing). The components of the 
vector <T and the order parameter are given by 
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where Ui are Pauli spin matrices, i — x,y,z, and A = 
AoiCTy. The above procedure is appropriate for describ- 
ing situations for which J <C Ep, w hich ho lds e.g. for 
the ferromagnetic alloys used in Refs. |l2ll,'^ . In writing 
Eq. we have chosen the normalization according to 



(4) 



Riccati parameterization: The spin-dependent na- 
ture of SF proximity systems calls for a formulation 
of the quasiclassical theory that retains the full spin- 
space structure, especially in studying situations where 
the exchange-field-orientation varies in space (such as 
in a domain wall). Within the Eilenberger theory a 
very convenient formulation already exists ^22 employing 
the so-called Riccati parameterizationiSiS The exten- 
sion of this method to the Usadel theory was achieved 
only recently, ■^'^ and has been applied to non-equilibrium 
situationSf^ and to FSF-systems with homogeneous 
magnetizations^ Here we demonstrate its usefulness by 



applying it to a SFS system with a spatially changing 
magnetization in a domain wall, a case where the con- 
ventional so-called ^-parameterization^ is not applica- 
ble. The spin-dependent Riccati parameterization)^ 
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automatically accounts for the normalization which 
is essential for practical numerical calculations. It is 
enough to determine one 2x2 matrix in spin space, 7. 
The other, 7, follows from the above-mentioned (funda- 
mental) symmetry. The transport equation for 7 follows 
from Eq. |5J, and reads^ 
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Here, the expression for J- is obtained by comparing 
Eq. O with Eqs. 

Boundary conditions: Additionally, boundary condi- 
tions are required for the different interfaces of the 
system. Such conditions have been formulated by 
Nazarovi^ The outer surfaces {z = Zo) of the supercon- 
ductors are assumed to border to an insulating region, 
and the appropriate condition is dzG(zo, E) — 0, i.e. 



^{zo,E) = 0. 
az 



(8) 



(3) On the other hand, the two inner SF interfaces {z 



for the S side, z 



for the F side) are assumed in 



the following transparent. The corresponding bound- 
ary conditions are in this case G{zf,E) = G{z[,E), 
csdzG{zf ,E) = (7FdzG{z[ , E), leading to 
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where as and ap refer to the conductivities of S and F, 
respectively. For simplicity, we have assumed 0-5 — up, 
implying the continuity of the derivative at the interface. 
With the boundary conditions ijSJl and lO, we have solved 
Eq. ITJ numerically by an iterative procedure (relaxation 
method) in the entire SFS system. 

SFS system with domain wall: We apply the outlined 
theory to study the SFS structure of Fig. in equilib- 
rium. Lengths are given in units of the superconducting 
coherence length, ^ = y/D/A^. The spin-singlet super- 
conductors are chosen to have the same gap magnitude. 
The contact areas at the SF interfaces are assumed to 
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FIG. 2: LDOS as a function of energy for the system of Fig.0 
calculated in the middle of the F region for different widths 
dw of the domain wall. Ifere, J — 0.5 Ao, dp ~ 2.0 For 
comparison, the dotted line shows the corresponding depen- 
dence in the normal-metal case (J = 0). The inset shows for 
J = 0.5 Ao, the value of the LDOS at the chemical potential 
{E = 0) as a function of the domain- wall width dw- 



be small enough, so that any gap suppression can be 
neglected. The two superconducting regions and the in- 
termediate ferromagnet are taken to have fixed lengths 
of ds = 5^ and dp = 2^. We model the domain wall 
by a varying direction of J = {Jx, Jy^ Jz) (keeping the 
magnitude J = |J| constant), with Jy = Q and 

Here, dw is an effective domain wall width parameter. In 
the following we study the influence of the width dw of 
a domain wall centered in F (zq = dp 12) on the density 
of states in the center of the domain wall [z = zo). 

Knowing 7(2, E) from the solution of Eq. Q with 
boundary conditions ||SJ| and Q , the quasiclassical Green 
fimction and the (total) LDOS 

Ntot{z, E) = --^ Im Tr g{z, E), (11) 

ZTT 

is determined via Eqs. (0)-®; Nq is the normal-state 
density of states, and Tr denotes the spin trace. 

An important characteristic of the value of the LDOS 
in superconductor-normal metal proximity systems is the 
minigap: the density of states in the normal-metal region 
shows a gap of width Eg < Aq induced by proximity to 
a superconductor. The energy Eg can be thought of as 
that of the lowest-energy Andreev bound state in a finite 
normal- metal layer. This convenient physical picture can 
easily be extended to single-domain ferromagnets. The 
corresponding spin-dependent energy shift of the quasi- 
particle and the Andreev-reflected quasihole by ± J leads 
to a reduction of the energy of the lowest-lying bound 
state, and correspondingly of the minigap, from the ex- 
pression for a normal metal by J, vanishing altogether 




FIG. 3: LDOS as a function of energy for several phase differ- 
ences (j) between the superconductors, calculated in the middle 
of F. Here, J = 0.5 Ao, dp = 2.0 The width of the domain 
wall is dw = 0.2 ^. The inset shows the corresponding LDOS 
at the chemical potential {E = 0) as a function of 0. 



when J > Eg j^Q. This picture is confirmed by our nu- 
merical calculations. 

In the inhomogeneously magnetized case of Fig.^ the 
above picture is modified. The effect of the domain wall 
on the LDOS is summarized in Fig. [3 which shows Ntot 
as a function of energy for different domain-wall widths 
dw- Although the value of J = O.SAq is here larger 
than the value of the normal state minigap Eg w 0.25Ao 
(as seen from the dotted curve in Fig. |21 for J — 0), 
the minigap is reduced to zero only for larger domain 
wall widths dw ~ 2^. For the smallest width dw = 
0.2^ the minigap is only reduced by about 40%. The 
additional states which fill the minigap with increasing 
domain wall width are due to spin triplet correlations, 
which are sensitive to the direction of J. Our calculations 
show that the influence of equal-spin pairing components 
created by the domain wall increases. This is reflected by 
the appearance of additional Andreev bound states inside 
the gap, modifying the LDOS. The relative importance of 
the triplet correlations depends on J and dw- as clearly 
manifested by Fig.|21 the efficiency of the triplet-inducing 
mechanism grows with increasing dw- The inset of Fig. [21 
shows the value of the LDOS at i? = as a function of 
dw- The interesting observation here is that the LDOS 
at the chemical potential is very sensitive to the domain 
wall width when the latter one is comparable to ^. 

Finally, with a view towards studying the possible ef- 
fects of the domain walls on the supercurrent flowing in 
an SFS structure, we have studied the LDOS in the case 
where there is a phase difference 4> between the two su- 
perconductors. This phase difference adds to the one 
accumulated by the quasiparticles and quasiholes in the 
ferromagnetic region and, thus, modifies the spectrum of 
Andreev bound states. Figures|31and0]present the LDOS 
in the middle of the F region for three domain walls with 
different widths. In the inset of Fig. |3|we also show the 
zero-energy LDOS for a domain wall of width dw = 0.2^ 
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FIG. 4: LDOS as a function of energy for several phase differ- 
ences 4> between the superconductors, calculated in the mid- 
dle of F, and for J = 0.5 Ao, dp = 2.0 The width of the 
domain wall is in (a) dw = 1.0 ^, and in (b) dw = 2.0 ^. 

as a function of the phase difference. As can be seen in 
Fig. ^ by a possible tuning of the domain wall width dw 
one can always find a region of strongest sensitivity for 
a given phase difference (f> and vice versa. This increases 
the possibilities of controlling the zero energy density of 
states in the domain wall. The rich structure exhibited 



by these results could easily result in highly nontrivial be- 
havior of the transport current both in equilibrium and 
nonequilibrium situations. 

Conclusions: We have studied numerically the LDOS 
in a heterostructure consisting of a ferromagnetic alloy 
sandwiched between two singlet superconductors. We 
find strong modifications of the LDOS caused by the 
presence of a domain wall. As only triplet supercon- 
ducting correlations are sensitive to the direction of the 
exchange field, the strong variations in the LDOS result 
from the presence of triplet correlations induced by the 
spatially varying magnetization. We also find a strong 
dependence of the density of states in the domain wall 
on a possible phase difference between the superconduct- 
ing order parameters, giving an additional tool to control 
its value. This motivates future studies of the interplay 
of a supercurrent and the domain wall (Josephson ef- 
fect). We hope that the variety of features observed in 
our calculations motivates further experimental research 
on proximity systems involving weak ferromagnets. 

This work was supported by Deutsche Forschungsge- 
meinschaft within the Center for Functional Nanostruc- 
tures. 
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